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Abstract 

We show that waveguides with a dielectric core and a lossy metamaterial cladding (metamaterial-dielectric 
guides) can support hybrid ordinary-surface modes previously only known for metal-dielectric waveguides. 
These hybrid modes are potentially useful for frequency filtering applications as sharp changes in field 
attenuation occur at tailorable frequencies. Our results also show that the surface modes of a metamaterial- 
dielectric waveguide with comparable electric and magnetic losses can be less lossy than the surface modes 
of an analogous metal-dielectric waveguide with electric losses alone. Through a characterization of both 
slab and cylindrical metamaterial-dielectric guides, we find that the surface modes of the cylindrical guides 
show promise as candidates for all-optical control of low-intensity pulses. 
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1. Introduction 

Electromagnetic metamaterials are materials designed to produce an electromagnetic response that is 
not possible with just the constituent materials alone [T]. Within the past decade metamaterials have been 
of interest [2l |3] , with metamaterials that operate at optical frequencies under current development [4H6] . 
Potential uses include cloaking devices [7], and perfect lenses [8]. Given that the electromagnetic response 
of a metamaterial is, in principle, tailorable, there has been considerable interest in developing metamaterial 
waveguides [5141 3| . The modes supported by metamaterial waveguides exhibit a wide range of properties with 
applications to improving the design of optical and near ultra-violet systems. We investigate the properties 
of modes in both slab and cylindrical waveguides with a lossy metamaterial cladding and show how both 
waveguide geometries are able to support low-loss surface modes. 

By exploiting effects that arise from surface plasmon-polaritons propagating along a single interface [21 
[T5] . a metamaterial-dielectric waveguide is able to support low- loss surface modes. The low- loss modes in 
a metamaterial-dielectric waveguide display a lower attenuation than is possible for the same mode in a 
metal-dielectric waveguide, thereby opening up the possibility of low-loss plasmonic waveguides. The low- 
loss surface modes we observe arise due to the same interference effects that cause the near-zero-loss surface 
mode observed on a single metamaterial-dielectric interface |14[ . The single-interface waveguide is used to 
confine the fields in order to enhance optical nonlinearity for all-optical control at a single-photon level [TS] . 
Rather than using a single-interface waveguide for all-optical control, employing a low-loss surface mode 
of a cylindrical waveguide, which confines the fields in the transverse direction, could enhance the optical 
nonlinearity that is required for all-optical control of pulses. 

Energy loss, which is a major obstacle in the implementation of metamaterial waveguides, is present 
in metamaterials due to scattering and heating effects in the materials used to construct them [THl. The 
artificial structure of metamaterials alters the scattering and heating effects of the materials [T7| . To fully 
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understand the properties of the modes in a metamaterial waveguide, it is crucial to include loss in the model. 
So far little consideration has been paid to the effects of loss [131 [18]. In order to model energy loss, we allow 
both the permeability and permittivity of the metamaterial to be complex-valued and frequency-dependent. 
For permeability we use a Lorentz-like model that describes the behavior of metamaterials designed for 
optical frequencies [17', whereas the standard Drude model is used for permittivity of the metamaterial. 

We compute dispersion, attenuation, and effective guide width for various modes of a cylindrical guide 
and, for comparison purposes, a slab guide. We show that a metamaterial-dielectric waveguide supports 
both surface and ordinary TM modes as well as hybrid ordinary-surface TM modes. Hybrid ordinary-surface 
modes are supported in metal-dielectric parallel-plate waveguides [19 . Hybrid modes are not predicted by 
lossless models. We find that a metamaterial-dielectric guide is able to support hybrid modes over a much 
wider frequency region than is seen in metal-dielectric guides. Some modes have attenuation curves with 
sharp changes over narrow frequency regions. The frequency at which this change occurs depends on the 
structural parameters, thereby enabling a metamaterial-dielectric guide to act as a frequency filter. 



2. Theory 

To characterize light in a waveguide, we solve the Helmholtz equation in each region (core and cladding) 
for both the electric field, E{x, y, z; t), and the magnetic field, H(x, y, z; t): 

V^F{x,y,z]Uj) = ^uj^e{Lj)^i{uj)F{x,y,z;uj), (1) 

where F{x, y, z; lj) is either E[x^ y, z; w) or H{x, y, z; w). We are interested in traveling plane- wave solutions, 
which have the form 

E{x,y,z;u)^Eix,y)e'^'''^-^'^ (2) 

and 

H{x,y,z;L,) = H{x,y)e'<-''''~'^'\ (3) 

for kz the complex propagation constant. As the permittivity and permeability of the claddings are complex 
functions, so too is the propagation constant kz := P + ia with a the attenuation per unit length for a 
propagating field and /? the propagation wavenumber. Both (3 and a are real and positive. A negative value 
for a implies the material is active (exhibits gain) whereas negative /3 corresponds to fields propagating in 
the negative z direction. Once the solutions to Eq. (jlj in each region are found, we apply the appropriate 
boundary conditions at each interface. Applying the boundary conditions yields the dispersion relation for 
the waveguide, which is geometry-specific. 

As we are interested in lossy waveguides, permittivity and permeability of the metamaterial must both be 
complex- valued and are denoted e = e' + ie" and fi = fi' + iiJ,", respectively, for e', e", and fi" all real, and 
e" > and fj," > 0. The expressions for e and /i are determined by how metamaterials are constructed |17] . 
A common approach used in metamaterial design for optical frequencies is to exploit the naturally occurring 
electric response of metals to achieve a negative permittivity and design the structure of the metamaterial 
to generate a magnetic response that gives a negative permeability [Sj. 

The natural electric response of metals to an electromagnetic field is described by a frequency-dependent 
expression for the permittivity given by the Drude model |20| : 

where We is the plasma frequency of the metal and Fg is a damping constant. The structurally induced 
magnetic response of a metamaterial takes the form (TT] 

^^1+ , (5) 
Mo uj^ - uj (uj + iT^} 

where F, a function of the geometry, is the magnetic oscillation strength, ujq is a binding frequency, and 
Fni is the damping constant for the magnetic interactions. The magnetic response is generated by resonant 
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Figure 1: A diagram of the structure of the slab and cyUndrical waveguides. Region (1) in both guides is the core and is a 
lossless dielectric, whereas regions (2), and region (3) of the slab guide are the cladding layers and are metamaterial. The width 
of the slab guide is w and the radius of the cylindrical guide is a = ui/2. Fields propagate in the z direction (into the page), 
whereas the x and y directions arc the transverse directions. 

currents that are driven by the electromagnetic field. For parameters that are real and positive, the real 
part of both Eqs. Q and ^ can be either positive or negative whereas the imaginary parts are strictly 
positive. 

The index of refraction of any lossy material is in general complex and given by 



where the sign is chosen such that ni > to ensure the material is passive (gain-free). The other case, ni < 0, 
corresponds to a material that exhibits gain; i.e., the intensity of the fields increases as the fields propagate 
through the material. The choice of the sign in Eq. ([6| to ensure rii > 0, means that at frequencies where 
both e' < and fi' < 0, the metamaterial behaves as a negative-index material with < 0. 

The frequency region where e' and /u' have opposite signs (one negative and one positive) correspond 
to metamaterials that do not allow propagating waves. Many metals at optical frequencies act as electric 
plasmas and are characterized by e' < and /x' > 0; hence, for metamaterials that are similarly characterized, 
we use the term metal-like to describe them. The metal- like region is what D'Aguanno et al. call the opaque 
region [T^- For most frequencies in the metal-like region rij > |n,.|. In contrast, metamaterials with e' > 
and /i' < are not metal-like per se but would correspond to materials that display electromagnetic responses 
consistent with an effective magnetic plasma [3]. For frequency regions with both e' > and fi' > 0, the 
metamaterial simply behaves as a dielectric with a frequency-dependent refractive index. 

We begin by studying the slab geometry as it has features in common with both a single interface (no 
transverse confinement) and cylindrical waveguides (supports multiple types of modes). Additionally, the 
slab waveguide has been studied for various materials [52] , which allows us to check our methods against 
accepted results. The slab waveguide consists of three layers as seen in Fig. [T] The center layer is the core, 
and the width of the waveguide, w, is the thickness of the core. We choose the width of the slab guide 
to be w = Attc/uJc, such that the frequency at which w is equal to the wavelength of the field is near the 
middle of the frequency range in which the metamaterial is metal-like. The core is surrounded on both sides 
by cladding. The cladding layers are sufficiently thick to substantially contain the evanescent fields, as the 
effective width of the guide (the width of the core plus the skin depth at the interface) is typically only 
slightly larger than w. The permittivity and permeability of the cladding material is described by Eqs. Q 
and (|5| whereas they are constant in the core. 

The slab guide supports two distinct classes of modes, transverse magnetic (TM) and transverse electric 
(TE). The TM modes have non-zero magnetic field components only in the transverse direction {H^^x = 0) 
whereas the TE modes have electric field components which are non-zero only in the transverse direction 
{Ez.x = 0). The existence of TM surface modes in our system is partially dependent on the cladding having 
e' < and /i' > 0, whereas the TE modes require e' > and /i' < in the cladding. As the expressions 




(6) 
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for e and /i, Eqs. Q and ([5| respectively, are different, we cannot expect TE surface modes to simply be a 
generalization of TM surface modes. 

We choose to study the TM modes as the forms of the response functions allow for the existence of TM 
surface modes for a large set of parameters, whereas the TE surface modes are not as robust. Applying the 
boundary conditions to the solutions of Eq. ([T]) for the TM modes leads to the dispersion relation 



ei / £2 £3 \ / ^1 , £263 
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tanh7ii/;, (7) 



which must be satisfied for the existence of a propagating mode . Here 



7j V fc? - ^^CjAij (8) 

are the complex wave numbers for the transverse components of the fields for j = 1, 2, 3 refering to the core 
and two cladding layers respectively. As Eq. ([7| is transcendental, we solve numerically for the complex 
propagation constant k^. To obtain the dispersion relation for the TE modes, simply replace any explicit 
appearance of €j in Eq. Q by /ij. 

Cylindrical guides enable transverse confinement of the modes but are of interest to us for more than this 
reason. The cylindrical geometry is well studied in other contexts and is used in a number of applications, 
e.g., fiber optics. Additionally, as for the slab geometry, the dispersion relation for a guide with cylindrical 
geometry can be obtained analytically. The cylindrical waveguide shown in Fig. [T] consists of a dielectric 
cylindrical core of circular cross-section with radius a. A cladding layer, which has a frequency dependence 
described by Eqs. Q and ([5]), surrounds the core. As with the slab geometry, the cladding layer is assumed 
to be arbitrarily thick. 

We again assume traveling-wave solutions, but we must include all components of both electric and 
magnetic fields. In general the modes of a cylindrical guide are not TE or TM but in-between modes that 
contain all of the field components and are called HE or EH [21]. The names HE and EH allude to the 
fact that all components of both electric and magnetic fields are non-zero. The TM and TE modes of the 
cylindrical guide are special cases of the general solution. Using the definitions 

J„ (a«;i) = — 

and 

Km (072) = 



dr 

where Jm and Km are the Bessel function and modified Bessel function, respectively, the dispersion relation 
for the cylindrical guide is [2T] 

1 , Y _ f fJ'l J'm i^Ki) ^ M2-^rn(Q72)^ f £l Jm {ai^l) ^ ^2 K'^{aj2) \ 



«1 I2J \K\Jm{aKi) ' Km{a'y2) J Jm{aKi) ' ^lKm{a^2), 

Here 72 is the complex wave number in the radial direction for the cladding and is given by Eq. ([s]) , and 
1^1 — ill is that for the core, whereas m is an integer characterizing the azimuthal symmetry. With the 
dispersion relations for both geometries, we are able to determine the propagation constants for the allowed 
modes. 



3. Modes 

The metamaterial-dielectric slab waveguide supports a number of TM modes with three distinct mode 
behaviors, namely ordinary, surface and hybrid modes. Each mode may display more than one behavior, 
with different behaviors supported for different frequency regions. The mode behavior is determined by the 
relative sizes of the real and imaginary parts of 71 , the wavenumber perpendicular to the interface inside 
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the core. To simplify the discussion we use ordinary wave, surface wave, and hybrid wave when referring to 
the character, or behavior, of a mode and use the term mode or exphcitly name the mode (e.g. TMj) when 
we are discussing a specific solution to the dispersion relation. 

We label the modes in a manner consistent with lossless waveguides [211: a numeric subscript j refers 
to ordinary waves and indicates how many zeros (how many times the Hy field changes sign) the mode has 
inside the core, whereas subscripts "s" and "a" indicate the symmetric and antisymmetric surface modes, 
respectively. Symmetric and anti-symmetric refers to the symmetry present in the Hy field. We extend the 
labelling scheme to the modes of lossy waveguides even though the defining features may not be present for 
all frequencies. 

Ordinary waves are caused by traveling electromagnetic fields inside the core reflecting from the core- 
cladding interface(s). The reflected fields overlap with the incident fields causing an interference effect, with 
constructive interference corresponding to the allowed modes. The condition for ordinary- wave behavior is 
Im (71) >> Re (71), which ensures that, in the slab guide, the fields have an oscillating wave pattern in the 
transverse direction of the core. 

Surface waves are the result of electrons at the core-cladding interface carrying energy along the waveg- 
uide. The energy from the fields is transferred to the electrons causing them to oscillate, providing a medium 
for energy transport. This type of collective oscillation of electrons at a surface coupled with electromagnetic 
waves is called a surface polariton or surface plasmon-polariton [23J . The polariton field of the surface mode 
decays exponentially as a function of the distance from the interface. Mathematically the condition for 
surface waves is Re (71) >> Im(7i), which causes the intensity to be concentrated around the interface(s) 
between the core and cladding and a relatively small intensity through most of the core. In slab guides only 
two surface wave TM modes are supported, namely symmetric (TMg) and anti-symmetric (TMa). 

In frequency regions where both Re (71) and Im (71) are comparable, the mode is a hybrid wave. Hybrid 
waves can be understood as a product of two features in the transverse direction, namely an evanescent 
feature and an oscillatory feature. To see how the two features combine, consider the components of the 
magnetic field in the core of the the slab waveguide, which have the form 

= iJo e'^i^ e^'^"^ e*('=^^-'^*\ (10) 

where Hq is a constant, the complex wavenumber 71 = 7^ + i7i , with ^[ and 7" both real, and the terms 
e^^^ and e*'''i^ represent surface-wave and ordinary-wave behavior, respectively. The hybrid waves are then 
a combination of surface waves, i.e., fields being carried by electron oscillations, and ordinary waves due to 
fields traveling in the core and reflecting off the interfaces thereby resulting in interference effects. 

For hybrid waves a substantial portion of the energy is transferred along the interfaces as well as in the 
core, so hybrid waves are like a combination of ordinary waves and surface waves. Hybrid waves are not seen 
in dielectric-dielectric guides and only for a very narrow frequency range of the TMa mode in metal-dielectric 
guides, where the mode changes from a surface wave to an ordinary wave. 

The metamaterial-dielectric cylindrical guide supports the same three mode types as the slab guide, 
though they differ slightly as the cylindrical guide has a circular symmetry. The fields of the ordinary wave 
TM modes in the cylindrical guide have the oscillating wave pattern of a vibrating circular membrane. The 
fields in the cladding decay exponentially, as Re (72) >> Im(72). The modes of the cylindrical guide are 
labeled similarly to the slab guide, with the transverse magnetic modes labeled as TM^ . For the cylindrical 
guide however, a numeric value of j refers to the number of zeros in Hy that lie between the center of the 
core and the core/cladding boundary. 

The HE modes have two indices (HE^ j), where the azimuthal symmetry is characterized by the integer 
m, which determines the number of oscillations in the flelds over an angle of 27r. The index m here is the 
same as that in Eq. The index j plays the same role for HE modes as it does for TM modes. Due 
to the circular symmetry of the cylindrical guide it supports only one surface wave TM mode. However, 
the cylindrical guide supports a number of surface wave HE modes of different orders with the symmetry 
determined, as with the ordinary wave TM modes, by the azimuthal parameter m. Modes of the cylindrical 
guide with surface-wave behavior are denoted by the index "s" for both TM and HE (e.g. TMg, HEi g. 
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Figure 2: Plots of the real parts (solid lines) and imaginary parts (dashed lines) of |(a)| e/eo (black) and ^//jq (gray), and |(b)| 
the refractive index of the metamaterial as a function of frequency. Regions (1) in each graph are the metal-like regions, region 
(2) is the negative index region and region (3) is the dielectric region. 



Re(yi)/Im(7i) ReCyO/ImCn) 




Figure 3: Plots of Re(7i) /Im(7i) for[(a)]the TMq (solid), TMi (dashed), TM2 (dotted), and TM3, (dot-dashed) and [(b)] 
the TMs (solid) and TMa (dashed) modes of the metamaterial-dielectric slab guide. The shaded region indicates the nominal 
range of hybrid waves, as there is no well-defined boundary. The insets are plots of \Hy\'^, which show the transverse profile of 
the mode at select frequencies. The thin vertical bars in the insets indicate the core/cladding boundary. 

4. Characterization 

With the physical nature of the three wave types supported by a metamaterial-dielectric guide under- 
stood, we now characterize the waveguides. To characterize the two guide geometries we examine how the 
effective refractive index, attenuation along the propagation direction, and effective guide width for a variety 
of supported modes vary with frequency. The following values are used for the metamaterial parameters |14] : 
Wo = 1-37 X lO^'^s"^ r„i = To = 2.73 x 10^^s"\ ojq = 0.2wc, F = 0.5 whereas the core is described by 
£1 = 1.3eo and /ii = ^q. 

Frequency regions with differing refractive index characteristics for the metamaterial are identifiable 
(Fig. [2]) when the above parameters are used along with Eqs. (|4| and ([S]). Above Wc the metamaterial 
behaves as a lossy dielectric with both e' > and fi' > 0. There are two metal-like regions, one for 
O.Suc < a; < cjo and one for uj < 0.2ajo- We say metal-like because the regions are characterized by e' < 
and /i' > 0, which is typical for metals at optical frequencies but not seen in dielectrics. In the frequency 
range 0.2^0 ^ to ^ 0.3a;o the metamaterial has a negative-index. The binding frequency, wq, in Eq. ([5| 
causes rapid changes in the refractive index (Fig. |^[b")] ) around oj — 0.2u}c- 

We characterize the slab geometry first, and, to provide a basis for comparison, we also include results 
for a metal-dielectric waveguide. We assume the slab guide is symmetric, meaning all of the parameters 
are identical for both cladding layers. The parameters we use to describe the metal are the same as for the 
metamaterial, with the exception that ^2.3 = Mo- To characterize the behavior of the modes we compare the 
real and imaginary parts of 71. A plot of Re (71) /Im (71) (see Fig. [s]) shows where the real and imaginary 
parts of 7i are of comparable magnitude, indicated by shaded regions, which correspond to hybrid waves. 
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Outside the shaded regions, either the real or imaginary part is dominant and the mode is a surface or 
ordinary wave, respectively. The insets of the plots also show the transverse profile of select modes at 
frequencies that lie in regions of different mode behavior. 

The transverse mode profile does not always present a clear distinction between ordinary, surface, and 
hybrid modes. As noted earlier, the fields of ordinary waves are generally distributed evenly across the core 
and have an oscillatory pattern, whereas surface modes are characterized by the fields being concentrated 
at the core/cladding boundary and decaying toward the center of the core. Hybrid waves, however, have 
characteristics of both wave types. One distinguishing feature of hybrid waves is a large amount of the field 
is present in the cladding layer (fat tails) . The two insets for hybrid waves in Fig. [S] clearly show the fat 
tails indicative of hybrid waves. 

Four modes are plotted in Fig. t^[a)| showing the regions of ordinary and hybrid-wave behavior for 
each. Hybrid-wave behavior is present for Re (71) /Im(7i) > 10~^, which translates to the TMq mode for 
0.2a;c <uj< 0.3wo, the TMi mode for 0.2wo < w < 0.45a;o, the TM2 mode for 0.2ujc < < 0.75wc and the 
TM3 mode for 0.2Ldc ^ 1^ ^ O.OSojc. The insets are plots of \Hy\'^ for the TM3 mode and show the difference 
in the transverse profile between the ordinary and hybrid waves. The inset for uj « 1.2wo is an ordinary 
wave (Re (71) /Im (71) « 10^"^), whereas the inset for uj « O.Swc is a hybrid wave (Re (71) /Im (71) « 10~^). 

Figure \'\[h)\ shows the TMg and TMa modes of the metamaterial-dielectric slab guide. Hybrid-wave 
behavior is present for 10~^ < Re (71) /Im(7i) < 50, which translates to the TMg mode for 0.2a;c ^ w < 
0.3wc and the TMa mode for O.lSwc ^ w ^ OAuJc- From the plots, we see that the TMg and TMa modes 
are surface waves for frequencies corresponding to Re (71) /Im (71) > 50. The insets in this figure are plots 
of \Hy\'^ for the TMg mode and show that when Re (71) /Im(7i) > 50 the mode is a surface wave (fields 
concentrated at core/cladding boundary), whereas it is a hybrid wave when Re (71) /Im (71) < 50. 

All guided modes in a waveguide experience an effective refractive index along the propagation direction. 
The effective refractive index of a waveguide is related to the propagation constant, /3, by the relation 
"off — P/ko with /cq = uj/c the free-space wave number, and c is the speed of light in vacuum. The effective 
refractive index of ordinary waves in a metamaterial-dielectric guide is very similar to that of a dielectric- 
dielectric guide, increasing with the frequency of the wave and bound from above by the refractive index of 



the core, rtcorc ~ jrv^elMi- The effective refractive index of the TMq and TMi modes is plotted in Fig. 4 [a.) 
For frequencies above 0.3cJe and 0.45cjo the TMq and TMi modes, respectively, are ordinary waves. 

Modes with surface-wave behavior exist in metal-dielectric guides as well as in metamaterial-dielectric 
guides. The requirement of e' or /i' to have opposite signs on either side of the interface means surface 
modes are not possible in dielectric-dielectric guides. Only two modes in the slab guide show surface-wave 
characteristics, TMg and TMa, which is present at frequencies such that ricff > '^coro- The effective refractive 
index of surface waves in a metamaterial-dielectric guide behaves as for a metal-dielectric guide, increasing 
with the frequency, and is plotted for the slab-guide modes TMg and TMa in Fig- ^'c) 



Hybrid waves sometimes exhibit anomalous dispersion (decreasing n^g with increasing frequency). We 
have already seen hybrid-mode behavior in Fig. ^|[a) where it is present at lower frequencies than the 



ordinary wave of the same mode. There are also regions of hybrid-wave behavior in the TMg and TMg 
modes, shown in Fig. ^[c) which occur for 0.2wo ^ w ^ O.Bojc for the TMg mode and for 0.2a;e ^ w < 0.4we 
for the T Ma mo de. 

Figure ^'b) shows the effective refractive index for the TM2 and TM3 modes of the metamaterial- 
dielectric slab waveguide. The hybrid-wave behavior occurs for 0.2wo ^ ^ 0.75a;o in the TM2 mode and 
0.2a;c ^ ^ 0.95ajo in the TM3 mode. Both modes have regions of anomalous dispersion with a particularly 
large frequency region of anomalous dispersion in the TM3 mode. 

The effective guide width of a waveguide, WcS, is the nominal width, w, plus the skin-depth at each 
interface, where = w + 1/Re(72) + 1/Re(73) for the slab guide. Ordinary and surface waves in a 
metamaterial-dielectric guide have an effective width that is only slightly larger than the nominal width; 
this is consistent with ordinary and surface waves of a metal guide. Hybrid waves, however, display an 
increased effective width when compared to the ordinary and surface waves. Figures f a) and |^[b)| are 



plots of the relative effective width Wcs/w and show the large effective width of the hybrid waves. The 
largest effective width for each mode is seen in the negative index region, where all modes are hybrid waves. 
Significantly increased effective width is a feature unique to metamaterial guides. The effective width of 
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metal-dielectric and dielectric-dielectric guides does not have such a large variation. 

A large effective width is due to the field having a large extent in the cladding, meaning a significant 
fraction of the energy is carried outside the core. As the cladding dissipates energy, there is a connection 
between effective width and losses. The surface waves, however, can have a small effective width even when 
there is a large fraction of the total energy in the cladding, which is due to the fact that the energy of a 
surface mode is concentrated around the interface with almost no energy in the center of the core. Thus, even 
when there is little penetration of the fields into the cladding, there can be a comparably small penetration 
away from the interface into the core. The effect, then, is a relatively even distribution of energy about the 
interface and a small effective width. 

The attenuation of a mode in a waveguide is a measure of how much the intensity of the mode decreases 
as a function of the distance travelled. The imaginary part of is the attenuation constant, denoted a, and 
characterizes the amount of attenuation per unit length that affects a propagating mode. All of the mode 
types in the mctamaterial guide show some amount of attenuation, with some wave types having more than 
others. 

Hybrid waves have the largest attenuation, which is partly due to their large effective width. The 
attenuation of the TM2 and TM3 modes of the slab guide is be seen in Fig. (^|[a) The regions of large 



attenuation occur where the modes display hybrid-wave behavior. Hybrid waves also have larger attenuation 
than surface waves in the same mode, though the difference is much less pronounced. The attenuation of 
the TMs and TMa modes is shown in Fig. |^[b)[ 

At frequencies where the mode displays hybrid-wave behavior, the attenuation is considerably larger than 
for either ordinary-wave or surface-wave behavior. The exception is the TMa mode for lu < 0.2ujc where it 
is an ordinary wave. Here the attenuation is large also but is consistent with the TMa mode in the metal 
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Figure 5: Effective guide width for | (a) [ the TMq to TM3, and |(b)| the TMs and TMa modes of the slab guide. Sohd lines are 
for the metamaterial-dielectric guide, and dashed lines are for the metal-dielectric guide. 
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Figure 6: Attenuation for |(a)| the TM2 and TM3, and |(b)| the TMs and TMa modes of the slab guide. Solid lines are for the 
metamaterial-dielectric guide, and dashed linos are for the metal-dielectric guide. 



guide (see the dashed hne in Fig. |f|b)), which is also an ordinary wave at low frequencies. The attenuation 
of both ordinary and hybrid waves decreases as the frequency increases, whereas the attenuation of surface 
waves generally increases with increasing frequency. 

As with the slab guide, we provide plots of Re (71) /Ini(7i) for the cylindrical metamaterial-dielectric 
guide in Fig.[7|to clarify where the modes are hybrid waves. The fat tails indicative of hybrid waves in a slab 
guide are also present for hybrid waves in a cylindrical guide. The transverse profiles of the hybrid waves in 
Figs. [3] and [7] show this increased field penetration, which is also reflected in the effective guide width for the 
modes. The plots in Fig. [7] show at which frequencies the modes of the cylindrical metamaterial-dielectric 
guide display which behavior. 

Figure 1 ^a) shows plots of Re (71) /Im(7i) for the modes HEio to HE12 and TMq. From these plots 
the different frequency regions exhibit where the modes display hybrid and ordinary-wave behavior. What 
is immediately clear is that there are no frequencies where the HEn mode is a hybrid wave. Instead it 
has a cut-off frequency near O.Swo. The three other modes (HEio, HE12 and TMg) are hybrid waves for 
< w < O.Swc- Beyond about O.Swo all four modes are ordinary waves. The insets in Fig. 7 ^a) are 

mode and show the transverse mode profile at two frequencies, uj sa O.Swe, where the mode is 
an ordinary wave, and lu ~ 0.25a;c, where the mode is a hybrid wave. 

Similarly, Fig. |'<||;b)| shows plots of Re (71) /Im (71) for the modes HE13 and TMi to TM3. Two of the 
modes (HE13 and TM2) do not show hybrid- wave behavior, whereas the TMi and TM3 are hybrid waves 
for Q.2uJe ^ w < 0.4aJe and Q.2u)c ^ w < 0.95a;o respectively. The insets here are plots of \H^\'^ for the TM3 
mode at two frequencies. Again, one of these plots represents an ordinary mode, lu w 1.2tJc, and the other 



for the HEio 



9 



Re(yi)/Ira(yi) 




Figure 7: Plots of Re (71) /Im(7i) for[(a)]the HEio (solid), HEn (daslied), HE12 (dotted), and TMq, (dot-dashed) , [(b)] the 
HE13 (solid), TMi (dashed), TM2 (dotted), and TM3, (dot-dashed) and the HEis (sohd) and TMg (dashed) modes of 
the cylindrical metamaterial-dielectric guide. The shaded regions indicate the nominal range of hybrid waves, as there is no 
well-defined boundary. The insets are plots of for half of the guide at select frequencies and show the transverse profile 

of the mode at these frequencies. The thin vertical bar in the insets indicates the core/cladding boundary. 



represents a hybrid mode, w ~ O.Swo- 

Modes with surface- wave behavior are supported by the metamaterial-dielectric cylindrical guide as well. 
Plots of Re (71) /Im(7i) for the HEis and TMg modes, both of which are surface waves for OAuc ^ ut < 
0.65ujc, are shown in Fig. I^c) These modes are hybrid waves otherwise, except for ui < 0.2ajo where they 
are ordinary waves. The insets in Fig. 7 ^c) show plots of \H^\'^ for the TMg mode showing a surface wave, 
Lu w 0.5ujc, and a hybrid wave, uj w 0.250;^. 

The modes of a cylindrical guide display similar dispersion characteristics to those in a slab guide, 
including the appearance of hybrid waves at lower frequencies. Figures I ^a) to (c) show nog for the modes 
HEio to HE13, TMo to TM3, as well as HEg, and TMg of the metamaterial cylindrical guide. In contrast to 
the slab guide, which supports only two surface modes, the cylindrical guide supports a number of modes 
with surface- wave behavior (one TM and the rest HE) , all of different orders (number of oscillations in the 
electric field over an angle of 2tt). 

The effective widths of the modes supported by the metamaterial-dielectric cylindrical guide show similar 
characteristics to the effective widths of the modes in a metamaterial-dielectric slab guide. Figure i |[a)| shows 
the effective width of the four modes HE13, TMi, TM2, and TM3 of the cylindrical guide, and Fig. [^[b)] 
shows the width for the HEis and TMg modes. The effective widths of the ordinary and surface waves are 
only slightly larger than the nominal width of the guide, whereas the effective widths of the hybrid waves 
are large, when compared to the nominal guide width and the effective widths of the other wave types. 

As with the effective widths, the attenuation of the modes of the metamaterial-dielectric cylindrical guide 
is similar to that of the slab guide. The attenuation for the HE13, TMi, TM2, and TM3 modes is shown 
In frequency regions of hybrid-wave behavior the modes have large attenuation. 



in Fig. IC a) In frequency regions of hybrid-wave behavior the modes have large attenuation. The same 
trend of decreasing attenuation with increasing frequency is also present. The TMg and HEig modes of the 
cylindrical guide, shown in Fig. l[ [b)[ have attenuation characteristics similar to the TM^ mode of the slab 
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guide. 

We now understand how metamaterial-dielectric waveguides behave when the effect of absorption is 
considered. The subject of possible imphcations and implementations of such waveguides can now be 
approached. 



5. Discussion: Low-loss surface modes 

In a metamaterial-dielectric guide, attenuation of ordinary and surface waves is much less than for hybrid 
waves but is still comparable to attenuation in metal-dielectric waveguides. With current metamaterial 
technology, losses due to atomic interaction with the magnetic field are significant |6l [17] . Previous results 
for optical metamaterials show magnetic losses could be reduced through structural improvements .24j . Our 
results show that lowering the magnetic losses opens up the possibility of low-loss metamaterial-dielectric 
waveguides. The attenuation of metamaterial-dielectric waveguides is reduced due to interference effects in 
the metamaterial [TH [T^ . 



The losses associated with the magnetic interaction are parameterized by T^. Figure 11 a) shows a 
comparison of the TMg modes for the metamaterial slab guide (Fm = Fe = 2.73 x lO^'^s""'^), the metal slab 
guide (Fni — 0), the metamaterial slab guide with reduced magnetic losses (F„i = Fe/3), as well as the 
metamaterial slab guide with greatly reduced magnetic losses (F^ = Fc/lOO). We choose Fm — Fc/3 to 
generate the plots as this should be practically achievable [53]. The attenuation dip seen in metamaterial- 
dielectric waveguides is related to that seen for surface waves on a single metamaterial-dielectric interface |14l 
I15j . The solution for a single interface, as well as the attenuation dip, is recoverable from the metamaterial 
guides in the limit of large guide width. 
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Figure 9: Effective guide width for |(a)| the HE13 mode, and TMi to TM3, and |(b)| thc flEis and TMs modes of the metamaterial- 
dielectric cyhndrical guide. The solid lines indicate HE modes and the dotted lines indicate TM modes. 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 




CjIcJi; 



Figure 10: Attenuation for | (a) [ the HE13 and TMi to TM3, and |(b)| the HEis and TMs modes of the metamaterial-dielectric 
cylindrical guide. The solid lines indicate HE modes and the dotted lines indicate TM modes. 



Close inspection of Fig. 11 '&) reveals that for uj > 0.45cl>c the attenuation of a metamaterial-dielectric 
guide with r,„ = Fc is slightly less than that of a metal-dielectric guide. Reducing F^ has the dual effect 
of further reducing the attenuation of the metamaterial-dielectric guide as well as broadening the frequency 
window in which the attenuation is less than that of a metal-dielectric guide. Similarly, the attenuation 
of the HEis mode of the metamaterial-dielectric cylindrical guide has lower attenuation than that for the 
metal-dielectric guide for uj > 0.45a;o. Reducing F,n for the cylindrical guide has the effect of reducing modal 
attenuation, as it does for the slab guide. Also, when F,n = Fc/lOO the attenuation of the the HEis mode 
of the metamaterial guide at w ~ 0.3we is 0.36 times that for the metal guide. 

Reducing the value of Fo, which parameterizes loss due to atomic interactions with the electric field, 
also has the effect of reducing the total attenuation of the TM and HE modes in a metamaterial-dielectric 
waveguide. However, reducing Fc for a metal-dielectric waveguide also reduces the attenuation of its TM 
and HE modes. Thus for a metamaterial-dielectric waveguide, reducing Fc alone does not yield any benefit 
over a similar reduction for a metal-dielectric waveguide. The benefit is in combining the reduced Fc with a 
reduced F„i to compound the two effects, which significantly reduces attenuation in metamaterial-dielectric 
waveguides. 

All-optical control of low-intensity pulses requires low attenuation and strong transverse field confine- 
ment for large cross-phase modulation. The cylindrical metamaterial-dielectric guide must meet the two 
requirements simultaneously to suffice for all-optical control. The first requirement, low attenuation, allows 
a low-intensity pulse to pass through the guide without being lost. By extending the results of Kamli, 
Moiseev and Sanders |14j . we have shown cylindrical metamaterial-dielectric waveguides can have lower 
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Figure 11: Attenuation for |(a)| the TMs mode of the slab guide and |(b)| the HEis mode of the cylindrical guide, and |(c)| 
the effective width of the HEis mode of the cylindrical guide. The solid line is for the metamaterial-dielectric guide, the 
dashed line is for the metal-dielectric guide the dot-dashed line is for a metamaterial-dielectric guide with Fm = Fc/lOO for the 
metamaterial and the dotted line is for a metamaterial-dielectric guide with Fm = Fe/3 for the metamaterial. 

attenuation than metal-dielectric waveguides when the magnetic losses of the metamaterial are reduced. 

The second requirement is strong confinement of the field in the transverse direction (normal to the 
core/cladding interface). Strong confinement allows the pulse to be confined to a smaller volume, thus 
increasing the local field intensity and enhancing the nonlinear response. As the low-loss modes of the 
cylindrical guide are due to field expulsion from the metamaterial cladding, the fields are almost entirely 
contained in the core. Thus, the metamaterial-dielectric cylindrical waveguide is a good candidate for 
enhancing cross-phase modulation between low-intensity pulses. 

Metamaterial-dielectric waveguides support modes with three distinct regimes, ordinary, surface, and 
hybrid ordinary-surface waves. The existence of the hybrid wave regime is a direct result of the waveguide 
dissipating energy. Models for metamaterial waveguides that do not dissipate energy (lossless) do not predict 
hybrid modes [HI [TT] due to the fact that the wavenumbers in the lossless model are allowed to be either real 
or imaginary, but not complex. Thus, the supported modes are either surface or ordinary waves, respectively. 
Modelling the waveguide with a dissipative material allows the wavenumbers to be complex, meaning the 
modes can have characteristics of both surface and ordinary waves. Hybrid waves in metamaterial-dielectric 
guides are much more robust than those in a metal-dielectric guide; i.e., they are supported for a much 
broader frequency range. The wider frequency range may make observing hybrid waves in practice easier 
for metamaterial guides than for metal guides. 

Modes with hybrid-wave behavior also have large attenuation relative to the other supported wave types. 
One possible use of hybrid modes is frequency filtering to remove unwanted frequencies from pulses. Pulses 
with frequencies outside of the hybrid wave region are allowed to propagate further, either as an ordinary 
wave or as a surface wave depending on the mode. Any pulse with frequencies that fall in the region of large 
attenuation propagate as a hybrid wave and are quickly extinguished. 
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6. Summary 

Our results show that reducing the energy loss associated with the magnetic interactions reduces the 
attenuation of the surface modes in mctamatcrial-dicilcctric waveguides. For r,„ < To the attenuation is 
below that for a metal-dielectric guide, at certain frequencies, despite the metal and metamaterials having 
identical expressions for the permittivity. Reducing Fm further reduces the attenuation in the metamaterial- 
dielcctric guide as well as broadening the frequency window of reduced attenuation. For example, with our 
choice of parameters, the attenuation of HEig mode of the metamaterial-dielectric guide can be reduced to 
as low as about 0.36 times that of the HEis mode in the metal cylindrical guide. 

With the use of structural improvements, or some other means of reducing Fm, it is possible, in principle, 
to construct metamaterial-dielectric waveguide devices with low losses compared to a metal-dielectric guide. 
With the capability of transverse confinement while supporting low loss modes, the metamaterial cylindrical 
guide is a good candidate for enhancing cross-phase modulation and ultimately all-optical control of low- 
intensity pulses. 

We have examined the modal properties of a slab and cylindrical metamaterial-dielectric waveguide using 
a model that includes energy loss. We have shown that the modes of metamaterial-dielectric waveguides 
support three distinct regimes, namely ordinary, surface and hybrid, with hybrid modes being a feature 
unique to models that include energy loss. Ordinary and surface modes are well known and have been 
studied previously. Hybrid waves, however, have not been previously examined for metamaterial-dielectric 
waveguides. 

Ordinary waves have the bulk of their energy distributed throughout the core and are the result of internal 

reflections and interference effects. Surface waves, on the other hand, have their energy concentrated at the 
interface(s) of the guide, which is due to the fact that the energy is being transported along the guide 
through the oscillations of surface electrons. Hybrid waves, as their name suggests, show a combination of 
the two mechanisms and thus their respective energy distributions. 

Hybrid waves typically have a large effective width when compared to the other mode types, which 
implies a large amount of energy is in the cladding. When the energy of a mode is carried in a cladding 
that dissipates energy, the result is attenuation of the field. Our model allows for a lossy cladding so the 
increased amount of energy in the cladding for hybrid modes is consistent with the increased attenuation. 
Within the negative-index frequency range of the metamaterial, all of the modes of the waveguide are hybrid 
waves and display the associated large attenuation and increased effective guide width. 

The sharp change seen in the attenuation curves could be advantageous as a frequency filter. A pulse 
propagating through a metamaterial guide would have the frequency components in the high-attenuation 
regions removed through energy dissipation, while the remainder of the pulse contimies to propagate. The 
frequencies at which the attenuation change occurs may be selected through the design of the metamaterial. 
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